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^ ; ABSTRACT 

■ We consider the interaction between radiation, matter and a magnetic field in a compact, relativistic jet. The 

entrained matter accelerates outward as the jet breaks out of a star or other confining medium. In some cir- 
cumstances, such as gamma-ray bursts (GRBs), the magnetization of the jet is greatly reduced by an advected 
radiation field while the jet is optically thick to scattering. Where magnetic flux surfaces diverge rapidly, a 
strong outward Lorentz force develops and radiation and matter begin to decouple. The increase in magne- 
tization is coupled to a rapid growth in Lorentz factor. We take two approaches to this problem. The first 
examines the flow outside the fast magnetosonic critical surface, and calculates the flow speed and the angular 
distribution of the radiation field over a range of scattering depths. The second considers the flow structure on 
both sides of the critical surface in the optically thin regime, using a relaxation method. In both approaches, 
we find how the terminal Lorentz factor, and radial profile of the outflow, depend on the radiation intensity and 
optical depth at breakout. The effect of bulk Compton scattering on the radiation spectrum is calculated by a 
\ Monte Carlo method, while neglecting the effects of internal dissipation. The peak of the scattered spectrum 

sits near the seed peak if radiation pressure dominates the acceleration, but is pushed to a higher frequency if 
the Lorentz force dominates, and especially if the seed photon cone is broadened by interaction with a slower 
£5 \ component of the outflow. 

Subject headings: MHD — plasmas — radiative transfer — scattering — gamma rays: stars 
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£> . Gamma-ray bursts involve collimated, relativistic outfl ows, as deduced from their rapid variabi lity, extreme apparent energies 
(which can exceed the binding energy of a neutron star: Kulkar ni et al.lll999t lAmati et afll2002l) . and the expected presence of 
\ non-relativistic m aterial surrounding the engine. Th e jet is heated as it works through this denser material, which may represent 
, a stel lar envelope dWooslevlll993l: |Paczvnski 1998), or neutron-rich debris from a binary neutron star merger (e.g. Dessart et al. 
7—1 ■ 120091) . As a result, a nearly blackbody radiation field may carry a significant fraction of the energy flux near the point of breakout, 
(f-) | We focus here on strongly magnetized jets that are driven outward by a combination of the Lorentz force, and the force of radi- 
, ation scattering off ionized matter . The acceleration of such a 'hot electromagnetic outflow' (|Thompsonll994tlMeszaros & Reesl 
CO 1 H997t iDrenkhahn & Spruit! 12001 iThompsonl [20061 iGiannios & Spruitll2007t IZhang & YarJBOl ll) . in which radiation pressure 
t— ( ' dominates matter pressure, has been treated quantitatively in lRusso & Thompson! (12012) (Paper I) in the approximation that the 
J" , poloidal magnetic field lines threading the outflow are radial and monopolar. The radiation field is self-collimating outside the 
. ' scattering photosphere, but may continue to interact with slower material that it entrained by the jet. In Paper I, the outflow was 
followed both inside and outside the fast magnetosonic critical point. The radiation force is especially important outside the fast 
point: even where the kinetic energy flux of the entrained charged particles is small compared with the magnetic Poynting flux, 
• they provide an efficient couple between magnetic field and radiation. The relative influence of the two stresses on the asymp- 
totic Lorentz factor depends on the radiation compactness. Generally, the importance of radiation pressure is enhanced by bulk 
relativistic motion at the photosphere. 

In this paper we generalize the calculation of Paper I to include non-spherical effects. A magnetized outflow experi- 
ences a st rong Lorentz forc e where poloidal flux surfaces in the jet diverge from each other f aster than in a monopolar 
geometry (ICamenzindl 119871: iLi et al.l fl992bT: iBegelman &XJ Il994t IVlahakis & Konigll I2003allb1: iBeskin & Nokhrind 120061: 
Tchekhovskov et al. 2009). In particular, a magnetized jet accelerates rapidly when it breaks out of the confining material 
(Tchekhovskoy et al. 2010). The simulations in that paper demonstrated the effect for a cold magnetohydrodynamic (MHD) 
outflow, but did not include the effects of radiation pressure and drag. 

The magnetization of a hot electromagnetic outflow remains modest inside its scattering photosphere, where the radiation is 
tied to the matter, and the radiation enthalpy contributes to the inertia. Our first task in this paper is, therefore, to examine how the 
radiation field begins to decouple from the matter when the jet material breaks out. We define a bulk frame in which the radiation 
force vanishes, by taking angular moments of the radiation field, and then track the proportions of the energy flux carried by 
matter, radiation, and magnetic field, at both large and small optical depths. 

Given the flow profile so obtained, the radiation spectrum is calculated by a Monte Carlo method. Here we focus on the effects 
of bulk Compton scattering, which provide a direct probe of the outflow dynamics. We neglect the effects of internal dissipation 
by various process such as MHD wave damping, magnetic reconnection, or shocks. 
The second principal goal of this paper is to obtain the longitudinal motion along a magnetic flux surface, taking into account 
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both the radiation force and the singularity in the flow equations which appears at the fast point. Our focus here is on the zone near 
and o utside the transparency su rface; previous efforts to calculate the effect of pressure gradient forces on relativistic outflows 
(e.g. IVlahakis & Ko nigl 2003b) have focused on t he optically thin regime. W e argued in Paper I that the effect of a magnetic 
pressure gradient driven by internal reconnection (Drenkhahn & Spruit 2002) has been overestimated, because it neglects the 
addition to the outflow inertia from particle heating and a strong non-radial magnetic field. 

Coupled wind equations for the matter Lorentz factor and angular momentum are derived in an arbitrary poloidal field geometry, 
restricted to the case of small angles near the rotational axis, but allowing for arbitrary relative flaring of the flux surfaces. The fast 
point generally sits close to the breakout surface of the jet. Our main simplification of the problem is to impose a particular shape 
for the poloidal flux surfaces, and not to solve self-consistently for the cross-field force balance. Two constraints are applied to 
the imposed magnetic field profile: that the rate of flaring is causal, and that the transverse component of the radiation force is at 
most a perturbation to the transverse Lorentz force. 

The plan of this paper is as follows. Section 2 reviews the acceleration of a relativistic MHD outflow driven by the differential 
flaring of magnetic flux surfaces, and considers the radiation transfer equation in the limit of small angles. Equations are derived 
for the acceleration of a steady MHD outflow outside its fast point, in combination with the radial evolution of the magnetization, 
radiation energy flux, scattering depth, and the frame in which the radiation force vanishes. These equations are solved in 
particular cases relevant to GRB jets. Section 3 presents a simple model of a spreading thin jet outside its photosphere, and 
derives the corresponding steady flow equations for arbitrary radiation force and magnetization. The effect of radiation pressure 
on the fast point is considered analytically, and numerical solutions for the flow both inside and outside the fast point are presented. 
Section 4 describes Monte Carlo calculations of the emerging radiation spectrum in both the causal jet model of Section 2, and 
the optically thin model of Section 3. Section 5 summarizes our results. The Appendix presents a derivation of the radiation force 
in a thin, transparent jet. 

2. FLARING, HOT MAGNETIZED JET: TRANSITION TO LOW OPTICAL DEPTH (MODEL I). 

We consider a stationary, axisymmetric outflow of perfectly conducting material that is tied to a very strong magnetic field. 
The outflow is also a strong source of radiation, which scatters off the advected light particles (electrons as well as positrons). 
Matter pressure gradients are neglected in comparison with inertial and Lorentz forces as well as the radiation force. 

We start by considering the exchange of energy between different components of the outflow. Deviations from radial motion 
are assumed to be small compared with the angular width of the photon beam: that is, the interaction between matter and radiation 




FIG. 1. — Geometry and approximate scale of the flow solutions for jet model I. 



is calculated assuming radial matter motion, but allowance is made for strong radial Lorentz forces driven by a small amount of 
magnetic field line flaring. The beam angular width is set, more or less, by the Lorentz factor of the outflow at its transparency 
surface. Here allowance is made for a finite optical depth of the magnetofluid. By taking angular moments of the radiation field, 
we track the difference between the Lorentz factor of the matter, and of the frame in which the radiation force vanishes, as the 
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matter is accelerated by a strong Lorentz force. This approach is suited to a single-component magnetofluid, but also allows for 
the presence of a second, slower component that scatters the radiation field into a broader cone, and plausibly is present in GRBs 
(Paper I). 

Given the complications introduced by a finite optical depth, we now consider only supermagnetosonic outflows. In a second 
approach (Section^, we account for non-radial matter motion and follow the flow across the fast critical surface, but restrict the 
calculation to low optical depth. The geometry of the model is shown in Figure [6] After being launched by the central engine 
(with angular frequency ft) the flow enters the jet zone along the rotation axis. We ignore the details of the acceleration while 
the jet is still very optically thick, and laterally confined. Our calculation begins a short distance inside breakout (at radius r*), 
by which point the flow is assumed to be supermagnetosonic. Outside breakout, transverse pressure support effectively vanishes 
and field lines begin to diverge differentially. The outward Lorentz force increases dramatically over a narrow range of radius, 
until a loss of causal contact across the jet forces the flow lines to straighten out, and the acceleration is cut off. Although the 
scattering photosphere could, in principal, sit anywhere in the outflow, breakout is associated with a large drop in optical depth. 
In our calculations, the photosphere therefore usually sits just outside breakout. Low optical depth at breakout does produce an 
interesting imprint of bulk Compton scattering on the emergent spectrum (Section|4|i. 

2.1. Exchange of Energy between Radiation and Magnetofluid 

We consider the flow along a poloidal magnetic field line 0/(r), starting at a large enough radius that the streamline sits well 
outside the light cylinder of the central engine. Deviations from radial motion are neglected, except in so far that they influence 
the radial Lorentz force. Then the outflow has a fixed total luminosity per sterad, including contributions from matter, magnetic 
field, and radiation, 

dL dEh dLp dL~> 

— = — - + — - + — Z= const. (1) 
dVL dVL dtt dtl 

Here 

1 dL k 2 „ Tc 2 dM 



r 



^=^-I>, = -- (2) 

is the kinetic energy flux of material of proper density p, poloidal (radial) speed v p , and Lorentz factor V. The poloidal Poynting 
flux is expressed in terms of the electric and magnetic vectors E, B by 

I dL P - E x B 

Substituting E = -v x B/c into the induction equation gives dB/dt = V x (v x B), where v is the fluid velocity. The steady 
solution to this equation involves the pattern angular velocity fl j of the magnetic field, which is constant along a poloidal flux 
surface. It relates the toroidal components of B and v via 

v<h — &lfrsin9f 

= -2 i J -B p . (4) 

v p 

Substituting this into © gives 

1 dL P . BpB^ 

-? =— HfrsmOf— — -. (5) 

r 2 dil 1 S 4tt ' 

Far outside the light cylinder, the outflow rotates slowly and the magnetic field is predominantly toroidal: <C v p ~ c and 
> \B P \. Hence 

^ ^rsinO f f£-. (6) 

It is useful to normalize all components of the energy flux to the poloidal mass flux, which is conserved along a poloidal flux 
surface in a steady MHD wind, dM/d$ p = Tpv p /B p = const. Assuming further that f>l, the magnetization becomes 



dLp/dil _ (n f r sin6 f) 2 B p B p 



(dM/dtyc 2 4ttc 3 Tpvp 



(7) 



The radius r* and the label * denote a position in the jet where the confining medium changes rapidly, e.g. the jet moves beyond 
the photosphere of a Wolf-Rayet star. (We will distinguish this from an inner boundary r,- for the jet integration, which typically 
is set just interior to the breakout radius.) Taking r„ ^> c/ftf, 

a _ (rsm9 f ) 2 B p ^ 

(Note that our definition of a differs by a factor V from the one used by Tche khovskov et ail 120091 in a similar derivation.) 
Defining the normalized photon luminosity by 

dL^/dQ 
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the equation of energy conservation ([T) can be written 

r-T*=-(o--cTt+n-n*). (io) 

Note that 1Z is related to the photon compactness and the scattering depth measured outward from radius r by 

*= - a ~(2-6)r 2 r es ^, (11) 
rmc J all 

where oj is the Thomson cross section and in is the material inertia per scattering charge. The coefficient here depends on the 
acceleration of the outflow, as can be seen from the expression for the optical depth of a (radially moving) photon 



V f$ 17 17 

s (r) = a es (r) / [l-/3(r')] — — - /3(r) = - ~ ^ 
' /3(r')r' c c 



(12) 



[see equation (1261 1 for notation]. The coefficient is ~ 2 when the Lorentz factor is constant and reaches ~ 6 for a linear growth, 
r oc r. 



i i i Mini i i 1 1 mi i i nil i i i in i i i iniij^ i i i nil 
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FIG. 2. — Radiation compactness at breakout for different ratios of Poynting to radiation flux: 1 (red), 100 (blue). Rescaled scattering optical depth r es *rj = 
10, 10 2 , 10 3 is plotted as solid, long-dashed and short-dashed respectively. The compactness drops below these curves beyond the breakout point. Strong radiative 
driving increases the Lorentz factor above the photospheric value T(r T ) if x* > T(r r ) 3 . 



2.2. Importance of Radiative Driving in Outflows with a Relativistically Moving Photosphere 

Jets of a high magnetization can encounter a scattering photosphere not too far outside breakout from a confining medium such 
as a Wolf-Rayet envelope or or neutron-rich debris cloud. Let us take a fiducial luminosity AirdL^/dQ, = 10 51 L51 erg s -1 and a 
deconfinement radius r* = 10 10 r* io cm. The corresponding compactness (fTTl i is \* = 1 x ^ & Ls\r~\ Q {m p /m). Since the Lorentz 
factor increases rapidly after breakout due to MHD stresses, we fix r and then consider the condition for a photosphere to emerge 
at a radius r T > r„. This corresponds to T 2 (r T )lZ(r T ) ~ lQ 7 (m p /m). For example, if the jet is hot and strongly magnetized, 
TZ ~ a ~ 10 5 , and pairs have largely annihilated within the bulk of the jet material, then the photosphere emerges at T(r T ) ~ 10. 

The radiation field is capable of driving a light baryonic outflow to a terminal Lorentz factor (see Section 2 of Paper I), 

roo~[r(r r ) X (r T )] 1/4 . (13) 

Moderately relativistic motion at the photosphere enhances and, as we now motivate, a larger radiation compactness. When 
the Lorentz force is taken into account self-consistently, can be greater or smaller than (|T3l , as we detail in this paper. 

An upper limit on the photon compactness is derived by demanding that the fluid be optically thin at breakout, T es (r*) < 1. 
Then the photosphere sits at r T < r», and the compactness (fTTl i at breakout is 

y» = — y(r T ) ~ 6— r z (r r ) ^ -. (14) 

A r* n (dM/dn)c 2 
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Hence Y^ oc T(r r ) 3 ' 4 in a jet of a fixed TZ. The radiation flux at breakout is below the Poynting flux if 

X* <6-r 2 (r r )a». (15) 

In Figure |2] we relate the compactness, magnetization and optical depth at breakout for different ratios of Poynting to radiation 
flux. 

The aforementioned hot jet moving at Y{r T ) ~ 10 at its photosphere, with a magnetization a ~ TZ ~ 10 5 , can be pushed by 
radiation pressure up to a terminal Lorentz factor ~(1 x 10 9 ) ~ 200. More relativistic material accelerated by the Lorentz 
force will, alternatively, feel a retarding force from the radiation field. 

2.3. Cold MHD Flow without Radiation Pressure 

To begin, we review the case where photons are absent, and assume a thin jet in which the magnetic field lines have poloidal 
angle 9j <C 1 . Only a small differential bending of the field lines is needed to push a cold magnetofluid to large T: their polar 
angle must deviate from conical geometry by S9f/9j^Y/a. A basic constraint on the rate of bending is obtained if the transverse 
component of the Lorentz force in the matter rest frame is limited to 



J'xB'-/3-(J'xB')| 

P 



< 



47rr/r' 



The prime ' denotes the rest frame, in which r/Y is a characteristic causal distance, and (3 = v/c. Then 



d(S9 f ) 



< 



1 



(16) 



(17) 



dr ~ T' 

y ~ 1 /T. This is seen in the dynamic cold MHD calculations of Tchekhovskoy et al. (2009), where strong 



so that typically 

Lorentz forces are generated in a narrow fan near the jet edge 

Conservation of magnetic flux implies that B p r 2 9fd9f = const. Hence, writing 9/ = 9f(r*) + 
becomes 

d(89 f /9 fif ) 
dln9f„ 

The change in the ratio of Poynting and mass fluxes can then be written 



1 da d 

a* dr d9f* 



K 



Y9 



9f = 0f* + 59 f, equation dgj 
(18) 

(19) 



The envelope function K(9f) ~ 1 away from the jet axis, and vanishes on the axis given the assumption of axisymmetry. 

2.4. Radiation Force 

Given the relativistic motion of the matter, the radiation field can be assumed to interact with it via Thomson scattering. In a 
frame where the matter moves with velocity (3c, and a photon has wave vector k = kk, a scattering charge feels a force 



pra 



rud . 



^ / (l-fB-t) k-f3Y 2 {l-f3-k 



dn. 



Here /(/i) = J dvl u is the spectral intensity integrated over frequency, and [i is the direction cosine [i = cos(#) = k- f. 
It is useful to define angular moments of /, 

F„ = 2tt / d^-fJ-Tm = 2tt / (A^)"/(Am), 



(20) 



(21) 



so that for a narrow beam, A/x ~ <^ 1, the radiation energy flux is approximately equal to Fq = TZYpc^. We may define a 
frame moving at Lorentz factor r eq (speed /3 eq ) in which the radiation field is nearly isotropic and the radiation force vanishes. 
Defining the bulk frame radiation energy density by u', one has 



/(A M ) : 



cu'/4ir 



[r eq (i-/3eq^)] 4 (i+2r2 q A^) 4 



(22) 



Substituting this into ( f2TT > yields the relations 



Fi = 



4Y 2 

^ x eq 



F ; 



Expanding the radiation force (l20l in A/i gives 



AYh 



[F -4r 4 F 2 ] 



Fi== T F 1 . 

eq 



o~tFo 



4Y 2 c 



Y 

- 1 - en 



(23) 



(24) 
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The main approximation here is that each field line experiences a small differential bending, so that the bending angle is small 
comp ared with (2A/1) 1 / 2 . This is consistent with rapid acceleration by the Lorentz force near the jet edge (Tchekhovsko v et all 
2010), e.g. 89 f ~ 6j — 0f <C 9j. In the context of GRBs, we can also assume that the flow has propagated far outside the light 
cylinder, so that <C 1 and the toroidal radiation force can be neglected. 

2.5. Transfer of a Narrow Photon Beam Near a Relativistic Photosphere 

We work with the transfer equation in the inertial frame into which the outflow is expanding; a prime denotes the matter rest 
frame. The radiation transfer equation is written (e.g. Mihalas 1978) 

= a es (S v -/„), (25) 

as 



where 



is the grey scattering opacity, and 



: r(i - /? M )c4 = ^ (i - $n) = Q ; s * (i - fto (26) 

fit p\V f i*) 



Su = ;S' = l - j / du'lL = l - T / du( l - Bu)l v (27) 

VXi-MP 2[r(i-/?/xW ^ v ZTii-WPJ 

is the source function in the isotropic scattering approximation. The Doppler relation between unscattered and scattered photon 
frequencies is v(l — /3ju) = v{l - (3fi), and 

l' v = [T{\-(iy)]%- rffl'^ri^ (28) 

[T(l-/3fj)] z 

are the usual transformations. Integrating over frequency gives 

^- = a es (S-I); S = Js v dv = -^^-^ J IQx){\- (29) 
Setting A/i — > 0, the path length is ds = dr//j, ~ ofr, and one has dA/i/dr ~ — 2A/i/r. Making use of 

JL F|I = -^±1> 1 1 ^(A M )"7+ 1 1 dfi(Afx) n f r , (30) 



gives 



and 



?i(^)=-4r^l 1 - 4r4 ^ F »=-4TW '-k F "' <31> 



1 d , a „ \ a es * / . T 



F = iW^2 1-rT K 02) 



r 4 ^ v 8TVA*) 2 \ r 4 q/ 

where we have set /3 — ?> 1 in the coefficient. These two equations, in combination with ( 1231 , allow us to evolve T eq near the 
scattering photosphere. 

The radial evolution of V and 1Z is obtained by differentiating ( [Tol l, substituting ( [T9T > and ( |3~Tl i, and expressing the radiation 
energy flux in terms of 1Z, 

dT - d (K\_dn dn _ a es « ( r 4 ^ 

r dr~ a *d6 f \T) K dr' dr ~ 4V\r/r.Y { T^J ' { ' 

2.6. Numerical Results 

Profiles are obtained for T(x), r eq (jc), R(x) and 6f(x) by integrating in the radial dimension. We have made the substitutions 
d/dOj — > (^gradient and K = 9f*/9j in (l33l l. This choice forces the field-line bending to zero near the center of the jet. The 
gradient angle <5# gra dient is of the order of 6j, but may be smaller near the edge of t he jet as it emerges f r om a confining medium. 
For example, the strong acceleration seen near the jet edge in the simulations of Tche khovskov et aTl (120101) is consistent with 
^gradient ~ 7 _1 ~ 0.16 f, it would presumably be reduced if the jet did not have a sharp edge. In our fiducial model we consider a 
field line anchored at #/* =0.1, and take the jet opening half-angle to be 9j = 0.2. To illustrate the effect of jet breakout on the 
flow parameters, we begin the integration just inside the breakout radius, r,- = 0.8r*. 

In our first set of integrations, the scattering optical depth is chosen to be large at the breakout radius. At this point, the radiation 
is still tied to the matter, r eq (r„) = T(r*), but the rapid MHD acceleration experienced by the flaring jet quickly forces a transition 
to low r es . The optical depth for a photon propagating radially from r = r* is given by equation (fT2l . and is thus unknown a 
priori. To impose a particular value of r es (r»), we first choose an approximate value of a es », evolve the equations of motion 
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FIG. 3. — Profiles of Lorentz factor V (blue), magnetization a (green), and radiation energy flux 1Z (red), in a flaring MHD outflow (cr* = 10 3 ) with radiative 
driving. The radiation force vanishes in a frame moving with Lorentz factor r cq (dotted blue). Scattering depth integrated from radius r to infinity (black curve) 
drops rapidly from r cs » = 10 as the flow accelerates. The gradient scale <5f9 gra( ]i cllt sets the relative degree of flaring between neighboring field lines. Left panel: 
Hi = 1, right panel: TZi = 1000. Other parameters: n = 0.8r„, T, = Y eqi = 10. 
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FIG. 4. — Same as Figure|3] but for t cs * = 1 (top plots) and r es » =0.1 (bottom plots). 
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and then iterate. The Lorentz force term in dT/dr is only valid outside the fast magnetosonic critical point, and so we take 

1 /3 

r(r,) > r c ~ oy . R adial integrations are done using a 5th-order Runge-Kutta algorithm with adaptive step size (see Sections 
7.3. 7.5 of lKiusalaasll20l0l) . 

Results are plotted in Figure [3] for an outflow with magnetization er, = er* = 1000, and both low and high radiation fluxes 
(IZj = 1, 1000) at the inner boundary. The action of the Lorentz force is concentrated at a small radius where the flaring is most 
severe, causing an increase in T that is initially much faster than linear. When the radiation energy flux is weak compared with 
the magnetic Poynting flux, the outflow experiences strong but logarithmic acceleration after breakout with V oc ln^ 2 x, a direct 
consequence of causally limited flaring, 58 f ~ 1/T [equation (f33tl . While the optical depth is large, the photon field is advected 
with the plasma, remaining nearly isotropic in the comoving frame. Once the optical depth falls below unity the photon field 
decouples and is free to self-collimate, so that r eq oc x. 

High radiation fluxes {1Z ~ a) force the flow back to the shallower profile T ~ r eq , even while the magnetic flaring grows 
stronger (1/T is larger). The acceleration zone is therefore widened in the radial direction compared with the radiation-free jet. 

Quite generally, we find that the terminal Lorentz factor is insensitive to the initial radiation energy flux. The growth in radiation 
energy flux outside the photosphere is therefore largely compensated by a further reduction in outflow magnetization. 

We also consider a low scattering depth at the breakout radius. In this case the matter is only weakly coupled to the radiation 
field, and the Lorentz factor is forced well above r eq a small distance outside the breakout radius (Figure 2). As in the case of 
higher optical depths, high radiation fluxes limit the growth of T and force it toward r eq . But the mismatch between V and r eq 
remains unless lZ>a. Energy is transferred from the magnetic field to the photons, resulting in a significantly modified spectrum 
(Section|4]i. 

Faster jet flaring, corresponding to a jet edge with 5f9 gla dient = 0.56*;, produces faster acceleration and terminal Lorentz factors 
closer to c*, but otherwise qualitatively similar behavior (Figure |5}- 

A novel effect becomes clear when the compactness is large: the magnetization can show a significant reduction, dropping 
significantly below 1Z and even V, and therefore resulting in a weakly magnetized outflow. This is caused by the strong radiation 
drag at small radius, which restricts the growth of T which allows for stronger jet flaring. 
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3. FLARING, HOT MAGNETIZED JET: TRANSPARENT FLOW 
ACROSS THE FAST CRITICAL SURFACE (MODEL II) 

The dynamics of the outflow can be calculated more precisely in the optically thin regime, where the radiation field is prescribed 
at an emitting surface (radius r s ) and passively collimates outside that surface. This allows us to study the critical point structure 
of the flow, at the price of neglecting the effects of multiple scattering. In the spherical case, the emission surface could, if one 
wanted, be identified with the physical surface of a star. But the model of a passively collimating photon field can also be applied 
to outflows that are already relativistic at the photosphere, including those with a jet geometry. In this section, we solve the wind 
equations for a steady, flaring jet which is optically thin but sub-magnetosonic at breakout. In contrast to the model presented in 
Section|2j here we prescribe the flaring profile in advance. 



Jet Zone Strong Acceleration Coasting 




10 7 cm ~ W n cm Logarithmic Radial 

Coordinate 



FIG. 6. — Geometry and approximate scale of the flow solutions for jet model II. 



The flow geometry is shown in Figure [6] The photon source radius sits in the confined portion of the jet, inside the breakout 
radius r*. After breakout, the optically thin flow is accelerated though the magnetosonic surface, whose location and shape are 
calculated self consistently (analytic approximations to the position of the critical surface can be found in Section l3~5l ). We follow 
the flow along a field line 6 fir), situated well outside the light cylinder, from just outside the source radius. 

Our solution to the flow below the breakout point formally is in the optically thin regime. Because the jet has already typically 
attained relativistic motion before breaking out, we can view the photon emission as arising from a virtual surface located below 
the physical photosphere, at a radius r,. e ff ~ r T /T(r T ). At a high radiation intensity, the matter is locked into the bulk frame 
defined by the photon field, T ~ T eq ~ (r — r s )/(9jr s ), around breakout. This means that the flow profile closely mimics an 
optically thick, radiation-dominated flow inside breakout, and we expect that our flow solutions should adequately represent the 
dynamics of a jet which encounters a photosphere at a radius r* ~ r T > r s . 

Our procedure is first to choose the poloidal field geometry and radiation profile, and then evolve the energy and angular 
momentum along the poloidal flux surfaces. A simple description of the photon field is possible when the jet geometry is locally 
spherical - that is, when the streamlines are conical inside breakout. This constrains the non-radial Lorentz force to vanish at 
r < r*, which in the small-angle limit can be written as 

i(JxB)^--\^=0 B.oce-K (34) 

c 4-nrO ad 

A jet with such a line current profile will de-collimate at r > r*, as the external pressure is removed. This decollimation leads to 
rapid outward acceleration of the cold matter entrained in the jet, even in the absence of radiative forcing. 

Oth er jet profiles are easily constructed and may be more natural: in the cold MHD jet calculation of iTchekhovskov et alj 
(2010), the confining surface has a parabolic structure inside breakout, transitioning to a conical geometry outside. Nonetheless, 
the rad ial Lorentz factor profiles that we obtain are (in the absence of radiation) very similar to those of ITchekhovskov et al.l 
d2010 ). and only depend on the magnitude of the differential flaring between magnetic flux surfaces. 
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We focus on the local dynamics within magnetic flux surfaces, taking into account the effect of radiation pressure. This 
longitudinal dynamics is sensitive to the relative flaring rates of neighboring flux surfaces, but not to the global profile of Poynting 
flux transverse to the jet axis. Although the transverse force balance is not explicitly taken into account, we do check that i) the 
degree of magnetic flaring is consistent with causal stresses; and ii) that the transverse component of the radiation force is weak 
compared with the transverse Lorentz force (so that the radiation flow is not strong enough to comb out the field lines into a 
conical geometry: see Appendix D). 

3.1. Jet Properties 

To construct an optically thin radiation field, we consider the simplest case of uniform intensity I = J l v dv at the emission 
radius r s , as we did in Paper I, but now restrict the sampling of the radiation field to polar angles 9 < 9j. The emission patch 
covers a small angular disk of area ir(8 jr s ) 2 (Figure[6]l, and the luminosity per sterad is dL 1:t /dft ~ TrOjl. This allows an analytic 
calculation of the radiation force acting on a particle of arbitrary Lorentz factor and direction, which is presented in Appendix A. 
This result generalizes the simpler angular moment formalism used in Section|2]and presented in equation (124-b . 

We showed in Paper I that if we normalize the photon intensity and the angular width of the photon beam by fixing i) the 
radiation force and ii) the relativistic frame r eq in which this force vanishes, then other quantities, such as the mean power radiated 
by an electron in its bulk frame, are nearly identical to those obtained from a radiation field that is isotropic at a relativistically 
moving photosphere. 

The radiation streams freely outward at r > r s , and its cone contracts with increasing radius. The size and orientation of this 
cone now vary with distance from the jet axis (in contrast with the case of a spherical emission surface; Paper I). There is generally 
a misalignment of the direction of peak radiation intensity with respect to both the radial direction, and the local flow direction. 
The alignment is strongest at a small but finite distance from the rotation axis, and produces a peak in the radiation force thereQ 

We normalize distances to r s , but measure the photon compactness (fTTt at the breakout radius, 

r Q f r s o~ T r., nOjl 

x=-; u = -!— x * = — r- —j- (35) 

r s c mc 3 (x»-l) 

In a GRB outflow, the photosphere generally lies outside the light cylinder of the rotating engine, so we take uj = VLjr s /c > 1 in 
our calculations. In this context, the magnetization can be most simply defined as3 

° = I i 3 (*/«!)• (36) 
Neglecting the radiation field, the energy and angular momentum per unit rest mass are given by 

and in a steady MHD outflow are conserved along field lines. They can be written in a dimensionless form, 

At = r G —^- ~ T + a; C = Tx0 f ^ f-^- = C m + C P . (38) 

uixOf B r u) l x6f B r 

3.2. Poloidal Field Configuration 

To incorporate a strong radial Lorentz force into the outflow, we choose the poloidal flux surfaces by fixing the function 
8f(x,6f*), where 6*/* is the polar angle at the breakout radius. The Lorentz force is large and positive if neighboring flux surfaces 
diverge from each other more rapidly than in a monopolar geometry. The effect of this differential expansion appears in the wind 
equations via the function 

At a \ dln9 f d ( d f*\ 2 a * 

We focus on the dynamics along a single flux surface, and so do not have to consider the angular dependence of B, •* = B, •(/*). The 
critical point structure of the longitudinal flow is insensitive to angu lar g radients in the flow magnetization. 
The flaring profile used in our calculations is described in Section l3~4l followed by the numerical results. 

3.3. Longitudinal Wind Equations 

We now consider the longitudinal evolution of the outflow variables along a magnetic flux surface. A radiation force (l20l is 
added to the Euler equation, which becomes 

pT\ ■ V(Tv) = — [(V • E)E + (V x B) x B] + ^F rad . (40) 
4tt m 

1 A tiny portion of the outer jet sits inside the light cylinder of the engine, and formally retains strong rotation of its field lines, which reduces the overlap 
between the radiation and fluid flows. This effect will, in practice, probably be eliminated by turbulence in the jet. 

2 This definition, following iMichel fl969h and lGoldreich & Julianl fT970h . differs in terms Ofy^/QfrsmOf) to cr= [c 1 dM / dQ)~ l dL P / dVL, the definition used 
in Section[2] and to a = (c 3 dM / dfiy 1 9^, where = Vj-rsinf/B^ is the advection rate of toroidal flux. We use equation J36t is independent of radius if the 
poloidal field is restricted to be purely radial, but the last two definitions are non-constant at Oiv^/QfrsinBf). 



The ^-component of equation d40i > is 
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Taking the dot product of equation (l40b with the unit poloidal field vector B p , defining the longitudinal derivative d/ = B p ■ V, and 
taking the small-angle limit, we have 

dl Tc 2 -^-d,C m = --r-^V 9 ' (r9 f B,) + JL- (41) 
B F mA 

Vpd,Cm = Tkf 1 + {r6f)J t- (42) 

Here C m is the specific matter angular momentum [equation (l3Ftl. and we have made use of the fact that the poloidal flow velocity 
\ p is aligned with R p . The Coulomb force only contributes to the transverse force balance and does not appear in equations (f4~TT > 
or (l42l i. Both of these features are easily derived by noting that the toroidal electric field vanishes in a steady, axisymmetric MHD 
wind (E ■ B = 0), which implies that \ p xB f = and E p ■ B p = 0. 

In Appendix A we calculate the radiation force (l20t in a thin jet, and express the poloidal and toroidal components in terms of 
dimensionless functions Rj, Pj, 

-2 -2 

Ff = X* —Rj(r, T, &); Ff = X * —Pj(r, T, /3 ). (43) 

Rotation of the photon field at the emission surface tends to reduce the azimuthal drag. It can be incorporated by modifying the 
/^-dependence of equations d43l . as is discussed in Appendix B, but is generally negligible when the outflow lies far outside the 
light cylinder (to ^ 1). 

A good approximation to the poloidal force can be obtained on field lines (xuS)~ l -C 0/ <S 0j, 

^^(^--sf} (44) 

in agreement with equation (|24j. The result for a spherical emission surface (Paper I) differs only in the absence of the factor $j. 
The vanishing of the radiation force occurs at a significantly higher Lorentz factor when the photon beam is collimated, 

r eq ^^ (45) 

up to a numerical factor of order unity as shown in Figure [TBI 

The deviation of the field lines from a purely radial direction is measured by A6b = Bg/B r , which we take to be small, so 
that B p = (1 + AO^y^Br ~ B r , v p ~ v r . As is detailed in Appendix C, the derivatives along field lines on the right hand side of 
equations ( |4TT > and ( f42l > can be written in the small-angle approximation. Ignoring the cross-field force balance then allows us 
to express (HTl i. ( f42l > as two ordinary differential equations, which can be re-written in terms of dT/dl and dC m /dl. The various 
terms on the right-hand side of these equations can be separated into purely magnetocentrifugal pieces (which do not depend on 
the radiation force), the direct radiation force, and a cross term: 



dr T' a +r'+r' dc m (O CT +0CW+CC)> 



dl flefS dl /Ltrff 



(46) 



where 



where 



and 



Asl+ ^£* ; „ s ( l+ m (l+A ^ + ?±^A (49) 

P r B r \ f J x B r A dl 



(J3 r Ty v (x9 f uj) 2 /3 r r \ PI I xuj0 f (jiT 




is the effective inertia. 
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3.4. Poloidal Field Profile 

We now prescribe the poloidal field profile outside the breakout surface, which, in a steady jet, also determines the poloidal 
streamlines. The profile inside r* is assumed to be straight and conical, Of* = 0f tS . A strong Lorentz force is obtained outside r* 
if the net change 0f iOO -Of* in polar angle is itself a growing function of Of*. A simple choice, that is asymptotically conical at a 
large radius, is 



0/0) = , e f* (, x * 
f * 56 \ x 



(x > x*). 



This connects smoothly with the inner cone if a > I. The net change in polar angle is determined b>0 SO, 



Vf,oo ~ Vf* 



The local change in the field line direction, relative to the total bend, is 

A0 B B e /B r 2ax* 



(51) 



(52) 



(53) 



"f~"f* "f~"f* X—X.jf 

The flux spreading factor works out to A = 2—0f*/6f for any field-line profile of the form (15TV In the absence of radiation, the 
Lorentz factor can be obtained by imposing energy conservation. It depends on the flaring profile of the jet via 



r~r*+<j*-(T = r*+<Tji-A- 1 ) = T*+o-*- 6f l 6f * 1 



lOf/Of* — 1 



; r, + <7* 



o f * I so 

\+20 f */50' 



(54) 



The acceleration tends to be more concentrated in radius for smaller values of the parameter a; in what follows a = 2. In Figure 
|7]we show sample field lines given by (T5TT > with strong flaring (SO = 0.3) for several values of Of*. 
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FIG. 7. — Sample field line profiles of the type 1511 with strong flaring (89 = 0.3) and breakout radius x* =r*/r s = 3 (dashed line). The emission radius bounds 
the inner grey zone. 
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3.5. Position of Fast Magnetosonic Surface 

When the flow speed surpasses the fast magnetosonic speed, radial magnetic disturbances are swept downstream and cannot in- 
teract with the part of the jet interior to the fast critical surface. The inertia of the electromagnetic field also becomes insignificant 
in the radial force balance, so that radiation pressure is relatively more important. We first consider how the position x c = r c /r s of 
the critical surface is modified by field-line flaring, and then consider the effects of radiatio n. The critical surface sit s at infinite 
radius only if the poloidal magnetic field is constrained to be radial and radiation is absent dGoldreich & Julianlll970l) . 

The critical surface is obtained by setting /i e ff = in equations d46l ). Retaining x* = 0, and assuming a 3> 1, this corresponds to 
r ~ cr 1 / 3 [a vanishing coefficient of dT/dl in equation (IC4H , and '5 = 0. When 50 f = O(0 f), the magnetofluid rapidly accelerates 

3 This is the analog of the parameter <50 Erll di ent appearing in jet Model I, used to approximate the derivative in equation 433 i . 
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outside radius r*, and so we can expand A ~ 1 near this radius: 

dA 6 f * f x*^- 1 ** 1 dA 

dlnl 59 V x J x x l 9fU) dl 

Here we have approximated (3$ ~ \/x6fU) <C 1. Then, for a = 2, 

(56) 

In the absence of magnetic-field flaring, the radiation stress forces the fast surface in from infinity (Paper I). Taking instead 
50 = oo but allowing for finite x*> the fast surface corresponds to T ~ c 1 / 3 and T^. ~ T' . Then 

x c 4cr 5 / 3 

— ~ (57) 

x* X *(x*-l) 2 uj 2 6 2 ft 

This differs from the spherical case (Paper I) mainly by the factor 1/0^*. At a very high compactness, the flow is tied to the 
collimating radiation field, and thus the critical surface is pulled in to where T eq (x,9f) ~ T c ~ er 1 / 3 . 
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FIG. 8. — Magnetosonic surface insi de a flaring jet at low radiation compactness (dotted blue, \ ~ ► 0) and high compactness [dotted red, x/cr <C 1 but still 
satisfying the bound given in equation ID3H . Left panels: Radius x c = r c j r s of the magnetosonic point as a function of field line footprint angle, for weak flaring 
(top) and strong flaring (bottom). Right panel: Two dimensional depiction of the magnetosonic surface (weak flaring). 

The fast surface is shown as a function of angle in Figure [8] for a breakout radius x* = 3. At low radiation compactness, this 
surface typically lies just outside breakout, x c > x*, in agreement with equation (T56b . As the compactness is increased, the critical 
surface can either move inward or outward, depending on the location where r eq = cr 1//3 . The critical surface is typically pulled 

inward near the rotation axis if the magnetic field is weakly flared, and also at larger polar angles if x* ~S> O/a 1 ^ 3 . Then its position 
follows equation d57| > until reaching the high-compactness limit at 



4a 4 / 3 

If alternatively the breakout radius is small, then the critical surface is pushed out by radiation drag, 



reaching its high-compactness limit at 




X* 



(58) 



(59) 



(60) 
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The deviation of x c toward large radius that is seen close to the rotation axis is due to a combination of effects: a reduction in 
the outward Lorentz force due to the weaker field-line flaring; and a mis-match between the radiation and matter flows driven by 
strong rotation. The first effect dominates at low x* ■ The change in critical radius at high \* can be estimated using equation 
( IA17b for r eq near the axis: 



rV3 



(%« w - 1 07 1 ) 



uj6 



f 



(61) 



3.6. Numerical Results 

We now examine the solutions to the wind equations (f46b-(l50l> that we have derived for a geometrically and optically thin jet. 
The singularity at the fast magnetosonic critical point, and the stiffness of the equations associated with large values of a and x< 
means that simple integration techniques such as Runge-Kutta are inadequate. To det ermine the position of the c ritical point and 
the flow solution inside it, we use the relaxation method described in Paper I (see also lLondon & Flannery| [T982). 
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FIG. 9. — Acceleration in a thin, strongly magnetized jet (<r» = 1000, u = 100, 0j = 0.2), along a field line anchored at 9f t =0.1. The radiation compactness 
at the breakout radius x* = r t /r s = 3 is varied: x* = 1, 10 3 , 10 5 , coiTesponding to the top to bottom curves on the left side. Left panel: 88 = 10 (weak flaring); 
right panel: 89 = 1; bottom panel: 89 = 0.3 (strong flaring). (The distance along the field line / differs little from the radial coordinate since the degree of flaring 
is small). 
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The inner boundary radius r, of the integration is chosen somewhat differently than in Paper I: we set it to twice the photon 
emission radius (x, = 2) because we only evaluate the radiation force where photons propagate at small angles with respect to the 
jet axis (requiring that xi — 1 3> Oj). The solutions for T and C„, obtained by an integration inside the critical point are required to 
be smooth near x,; avoiding sharp gradients restricts the boundary values at x; to a narrow range. We also make a first guess for 
the critical point radius x c . The regularity of the solution at x c then allows us to determine the flow variables at the critical point 
from the equations 



fr' +r' +V) 

\ c °X X) 



= /leffO c ). 



(62) 



An approximate solution is chosen which connects the inner boundary values to the critical pointQ This solution, along with 
the position of the critical point, is then relaxed to within a desired tolerance using a Newton-Raphson method, all the while 
satisfying the regularity condition (l62l . As a last step, the flow outside the critical point is obtained by shooting outward using a 
fifth-order Runge-Kutta algorithm. 

Solutions are obtained for a range of photon compactness and a high magnetization (er* = 1000). The part of the jet studied 
sits well outside the light cylinder, u = 100 in equation (l35l ). Choosing the flaring profile (IBTt . we follow the flow along a 
field line with initial footprint ffy* = 0.1, in a jet of half-opening angle Oj = 0.2 and a breakout radius x* = 3 = 1.5x,. The 
magnitude of the jet flaring is adjusted by choosing the parameter SO, with values 10, 1,0.3 corresponding to a net angular shift 
Of, co I Of* - 1 = 0.10/*, 0/*, 3.3 Of ^ between breakout and infinity. The maximal flaring chosen (50 = 0.3) still satisfies equation 
( [17] ). and so the divergence of neighboring magnetic field lines is consistent with causal stresses. 

The results are show in Figure |9] At low radiation compactness, they resemble those obtained by Tchek hovskoy et aTl (120101) 
for a cold MHD jet. A slow, nearly linear, increase in Y within the star is followed by rapid (but logarithmic) growth beyond 
the breakout point, where the field lines begin to diverge. As \* increases above ~ a, photon drag begins to dominate the weak 
Lorentz force inside the breakout radius, and Y tends to T eq ~ x/Oj [equation d45l )l. After breaking out, the fluid is quickly 

1/3 

accelerated through the fast point. A strong radiation field forces the position of the critical point to a radius where r eq ~ oy (in 
this case, the displacement is outward). We do not search for solutions with the fast point inside the star, corresponding to x* > 3. 
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FIG. 10. — Asymptotic Lorentz factor of flaring magnetized jet (cr* = 1000, u = 100, 9j = 0.2), as a function of compactness x* at me breakout radius. Curves 
correspond to a field line anchored at 8f* =0.1 and represent varying degrees of flaring (weak: 69 = 10; intermediate: SB = 1; strong: S0 = 0.3). Dashed black 
line: asymptotic Lorentz factor 0.8(x*^*/46y) 1//4 of an unmagnetized, radial, monopolar outflow. 



The influence of the radiation field can be clearly seen in the dependence of asymptotic Lorentz factor Too on \* (Figure ITOb. 
When flaring is strong (50 < 1), the matter is pushed rapidly to Y > r eq and it feels a net drag outside the magnetosonic point. 
The asymptotic Lorentz factor is reduced from (|54V At very high x*> the radiation drag is able to suppress the acceleration and, 
therefore, the asymptotic Lorentz factor. 

The minimal compactness needed to significantly affect the post-breakout flow can be estimated by equating the leading terms 

4 This procedure is followed for a discrete set of values of x*> a ft er which the inner boundary condition on T, C m is chosen from an interpolating function. 
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in ( PTTl i at x = 3x*/2, the point where the relative flaring of poloidal field lines is maximal. This gives 

5/3%* 

X*>°* -jg- (63) 

In this case, the sign of the effect of radiation pressure outside the critical point still depends on the degree of magnetic field line 
flaring. If flaring is weaker (80 > 1) then the details of the flow profile near breakout differ from the unmagnetized flow, but 
is still well approximated by an unmagnetized, radiatively driven flow: 

Too -O.S(x*x*/0j) l/4 (cr = 0, 69 -> oo). (64) 

Our solutions, in the part of parameter space that we have explored, satisfy two basic constraints. First, the outflow is optically 
thin at (or near) breakout if the compactness sits below the bound ( TT3T >. Second, the component of the radiation force transverse 
to the poloidal flux surfaces must remain sma ll compared with the transverse Lorentz force that is implied by the chosen flaring 
profile. The corresponding upper bound (ID31 on the compactness is derived in Appendix iDl 

4. SPECTRUM OF SCATTERED PHOTONS 

We now consider the self-consistent spectrum of photons that scatter off a hot electromagnetic outflow near its photosphere. 
Our focus is solely on the signature of the differential flow of matter and photons - that is, we neglect any internal processes that 
would heat particles or induce small-scale deviations from a uniform flow. 

We first consider the jet model of Section [2] in which the optical depth is finite but t he fl ow is considered only outside its 
fast magnetosonic surface. Then we turn to the optically thin jet model of Sections [31 and 13.61 in which the entire flow is solved 
inside and outside the critical surface, but the region interior to the photosphere is ignored. As in Paper I, we neglect any 
internal dissipation in the outflow, which can contribute to the high-energy tail of the spectrum (Thompson 1994; Gianniosl l2006l : 
lBeloborodovl2010l) . 

4.1. Spectrum in Jet Model I ( Super-magnetosonic): Monte Carlo Method 

Here we follow the photon field self-consistently across the jet photosphere, which is assumed to sit outside the fast magne- 
tosonic point (FigureHJ. The exchange of energy between photons and magnetic field was calculated in Section[2]in parallel with 
the flaring rate of the poloidal field lines outside a fixed breakout radius. The radiation force on the matter, and the evolution of 
the equilibrium Lorentz factor r eq of the radiation field, defined in equations (l23l and d24l >. are both calculated by taking angular 
moments ( f2TT > of the intensity. 

To calculate the emergent spectrum, we i) take the flow velocity profile as a given background, and then ii) inject photons from 
the inner radius r, = 0.8r» with an isotropic distribution in a frame moving with Lorentz factor r eq ,,-. Photon parameters in the 

rest frame of the 'star' at r > r, are obtained by a simple Lorentz boost. Defining a radial direction cosine by p = k - r, one has 
[i = (// + /3)/(l+/5//), u> = T(l + f3p!)bj' , where the prime labels the matter rest frame. Deviations from radial flow are assumed 
small compared with the width of the photon beam. 

Electron scatterings are handled in the Thomson approximation (the outflow moves relativistically in the case of a GRB), by 
drawing a random number 1 -e~ Arc \ The position of the next scatter point is calculated by integrating 

[ n r*dr 
AT es [ri,r 2 ,/x(ri)] = a es *r* / [1- p(r)/3(r)] — (65) 

along the photon ray. Note that the outflow solution of Section [2] has been iterated so that the coefficient a es * corresponds to a 
prescribed value of the radial optical depth T es (r» , oo, 1) at the breakout radius. Other initial flow parameters are defined at r,-. The 
direction cosine evolves from a scattering radius r to r 2 > r according to 
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1-Aifor=(-J (r 2 >r). (66) 

The photon escapes if Ar es exceeds the total optical depth along the ray. 

The frequency distribution of the outgoing photons is first obtained with a monochromatic photon source, l u = Iqvq5(v- vq). 
This output spectrum is then convolved a source spectrum that is either a pure blackbody, or a function^ that mimics the low- 
frequency slope of a GRB, F v = const x e~ hu l kT °. For both types of seed spectrum, the temperature 7b is normalized by requiring 
F v to peak at the fixed reference frequency vq. 

Scatterings are taken to be elastic in the bulk frame, where the matter is assumed cold, so that the outgoing and ingoing 
frequencies satisfy the usual Doppler relation, 



V 1 -/J /iem 1+/3A*' 

After transforming p to the local bulk frame, we pick scattering angles 9' s , (f)' s with respect to the flow direction. The direction 
cosine of the outgoing photon is determined via p' em = /j,' cos8' s + ( l -/i' 2 ) 1 / 2 sin 0' s cos <j/ s , followed by a boost to the stellar frame. 



5 This corresponds to the lower-frequency half of the Band function I Band et al. 1993), extended to all frequencies. 
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The peak of the seed photon distribution is stretched to higher frequencies when the outflow Lorentz factor T > r eq [equation 
(|45j]. A scattered photon has a frequency in the range < v < v mdLX , where 

*W = 1 vq ^ U + ^ (x,T»l) (68) 

and vnfr = [(1 -/?)/(! +/3)]po ~ vo/4T 2 . 



4.2. Spectrum in Jet Model I: Results 
Figure[TT]shows spectra for the case where matter and radiation field are initially locked together at T, = r eq ,, = 10. The curves 
correspond to a variety of optical depths, as well as low and high initial photon fluxes, 1Z, ■ = 1, 10 3 . The peak of the spectrum 
is somewhat broadened compared with a pure blackbody, and the segment shortward in frequency of the peak has a flattened 
spectrum, although not as flat as is seen in GRBs. A similar effect was seen in Paper I in the case of hot electromagnetic outflows 
accelerating along a radial, monopolar monopolar magnetic field. The spectrum below the peak is flattened even more if the 
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FIG. 1 1 . — Spectra of photons emerging from a relativistic jet with a radial profile calculated by the method of Section|2] for various values of the optical depth 
at breakout t cs * = 0.1, 1, 10 [equation 1121 1. Black curves correspond to a monochromatic seed, green curves to a blackbody seed, red curves to a GRB-like seed 
spectrum, F„ = const X e~ hv l kT ° [the lower-frequency half of the Band (1993) function extended to all frequencies]. Normalized photon energy flux at the inner 
radius r; = 0.8r*: TZi = 1 (left), IZi = 10 3 (right). Initial Lorentz factor T, = 10. Bulk frame of the radiation field r K] ,, = 10, except for the dashed curves which 
correspond to r c q., = T//2 = 5, t cs * = 0.1. 
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radiation field emerging at the jet photosphere is broader than the matter Lorentz cone: the dashed curve in the IZj = 1 panels is 
the result for r eq ,,- = r,/2 = 5 and a low optical depth (r es » = 0.1) at the breakout radius . As expected from the above argument, 
the peak of the spectrum is stretched upward in frequency above the peak of the seed spectrum. 

4.3. Spectrum in Jet Model II (Trans-magnetosonic) 

The output spectrum is calculated by a simila r me thod to that described in Section 14.11 The background flow is prescribed, 
in this case by the solutions obtained in Section [3~6l and the flow is approximated as radial. Since we are not, now, following 
the outflow across its photosphere, and the entire simulation volume is assumed optically thin, we take a similar input photon 
distribution as was used to calculate the flow acceleration: the intensity is constant, / = /q for 9 <9j = 0.2 radian. The jet breakout 
radius is taken to be x* = 3, as in Figure|9] and the radial optical depth Ar res (x; = 2,oo,/i = 1) = 1. Results are shown in Figures 
IT2HT31 for three different degrees of magnetic field flaring: strong [corresponding to 59 = 0.3 in equation (BTTll. intermediate 
(59 = 1), and weak (59 = 10). The high-energy extension of the spectrum becomes broader as the radiation compactness is 
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FIG. 12. — Photon spectram emerging from a highly magnetized jet with radial profile given in Figure|9] corresponding to strong magnetic flaring [68 = 0.3 in 
equation 1511 1. The radiation compactness x* = 1. 10 3 , 10 at the breakout radius x* = 3. (The spectra extend to higher frequency at a lower compactness.) The 
optical depth to radially moving photons is unity at the inner boundary x = 2. Black lines: monochromatic seed spectrum. Left Panel: black-body photon source 
(green lines). Right Panel: GRB-like seed spectram, F v = const X e~ hv ' kTa . Dotted curves: source spectram. 



reduced: stronger radiation drag limits the increase in T above r eq . In the case of a monochromatic input spectrum, one notices 
the appearance of a few distinct orders of Compton scattering. The bumps in the spectrum are smoothed out when convolved 
with a blackbody source. 

5. SUMMARY 

We have examined the effect of intense radiation pressure on a cold, magnetized outflow with a jet geometry. The poloidal 
magnetic field lines are allowed to deviate from spherical symmetry, e.g. due to breakout from a confining medium. The outflow 
experiences a strong outward Lorentz force as a result, so that the magnetofluid and radiation field have a tendency to flow 
differentially outside the transparency surface. We have considered the combined dynamics of the magnetofluid and radiation, 
and as well as the modification to the radiation spectrum by multiple scattering. 

We first considered the transition zone straddling the scattering photosphere. While the jet is still optically thick, its magne- 
tization is suppressed by the inertia of the advected radiation. Outside the breakout point, the jet experiences a strong outward 
Lorentz force, which forces a rapid reduction in optical depth. This approach assumes that the fast critical surface lies deep in 
the jet, but calculates the radial flaring of the field lines self-consistently with a simple causal prescription, and calculates the 
interaction between the radiation and matter for arbitrary scattering depth. 

If the jet is still optically thick at breakout, then the emergent spectrum is modestly broadened and hardened below the peak. On 
the other hand, breakout outside the transparency surface results in a photon beam that is significantly broader than the Lorentz 
cone of the accelerating jet, and therefore results in a more extended high-energy component to the spectrum. A stronger radiation 
field suppresses the accelerating effect of jet flaring and brings the spectrum closer to the original the rmal input. Broad ening of 
the photon beam could also be due to scattering by a shell of slower material entrained at the jet head (lThompsonl2006l) . 

Our second approach to the problem focuses on the zone outside the jet photosphere, but allows for large enough magnetization 
that the flow passes through the fast critical surface just outside the breakout radius. We then solve for the flow profile along 
magnetic flux surfaces, both inside and outside the critical point. In doing this, we choose a realistic angular distribution for 
the radiation field but prescribe a flaring profile for the poloidal field lines. The cross-field force balance is not solved self- 
consistently, but we check that in all cases the transverse component of the radiation force is small compared with the transverse 
Lorentz force that is implied by the chosen field profiles. 
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FIG. 13.— Same asff2] but with 56 = 1= 56 j (top), and &6 = 10 = 506*/ (bottom). 



As regards the longitudinal motion along magnetic flux surfaces, we define the critical compactness \ of the radiation field 
above which the matter and radiation are locked, and the Lorentz force is subdominant. For small jet flaring, the radiation force 
leads to an increase in terminal Lorentz factor at high values of %, but can somewhat suppress the acceleration if the flaring is 
strong. The extent of the high-energy component of the spectrum is shown to depend in an interesting way on the degree of 
flaring and the position of the photosphere relative to the breakout radius. 

Issues not addressed in this paper include the eff ects of multiple scattering at the m agnetosonic critical surface, an ambient 
radiation field generated far outside the engine (e.g. iLi et al.| [l992a; Bes kin et alj|2004l) . or the feedback of an intense radiation 
flow on the poloidal structure of the magnetic field. An effect specific to gamma-ray bursts involves the sidescattering of gamma- 
rays outside the forward shock, combined w ith the radiative acceleration of the pair-enrich ed material up to a Lorentz factor 
comparable to that of the relativistic ejecta (Thompson & Madau 2000; Beloborodov 2002). This delays the deceleration of 
the ejecta, and makes the medium ahead of the shock optically thick to scattering (lThompsoni r2006). Photons side-scattered 
through large angles would continue to interact with jet material at a smaller radius, creating pairs downstream of the forward 
shock, delaying the decoupling of the photons from the jet fluid, and generating a high-energy tail to the photon spectrum by bulk 
Comptonization. This means that the outermost shell of jet material (of a thickness ~ &] r *) ma y avoid strong outward acceleration 
during jet breakout. However, jet material flowing at much greater distances back of the jet head sees weaker Compton drag and 
rates of pair creation during breako ut. The slow forw ard shell becomes geometrically thin as it is pushed outward and, eventually, 
subject to a corrugation instability (Thompson 2006). 
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APPENDIX 

GEOMETRY OF SCATTERING IN A NARROW JET 

We now calculate the radiat ion force on plasma moving on a general trajectory within a thin jet of opening angle 0j <C 1, 
following the setup of Section [XT] Our goal is to obtain an analytic expression for this force, which is possible by assuming a 
uniform intensity 1 = J l v dv at the 'emission' surface (radius r„), and taking this surface to be locally spherical. When considering 
the interaction with matter, this intensity distribution gives similar res ults t o a radiation field that is locally isotropic in the 
relativistic frame of the emitting medium. (See the discussion in Section [37X1 and Section 3.3 of Paper I.) The result generalizes 
the simpler angular moment formalism used in Section |2]and presented in equation ( f24b . 

Photons are emitted from coordinates {# 7 ,(?!> 7 } within a patch of angular radius 1 < Oj, and scatter in the jet the position 
[x = r/r s > 1, Of, 4>f = 0}. The presence of an absorbing surface at the edge of the jet would change the radiation force at angles 
Of > Oj. Given the uncertain nature of the medium outside the jet, we restrict the calculation of the force to angles < Oj. The 
photon trajectory is tilted with respect to the radial line passing through the scattering point, by an angle (Figure \14] 

7 , r =^f. (Al) 
x- 1 

Here 

(Q ir ) s = (e 2 f + 1 1 -20fe i con(j) 1 ) 112 (A2) 

is the corresponding angle measured on the 'emission' surface, and we have assumed that x— 1 ^ Oj in making the expansion in 
x. The intensity at the point of scattering can then be expressed as 



/= / dvl v =h for 1 . r < Jo 2 f + 2 j -20f0jcos(j) 1 . (A3) 
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The unit wave vector in the local (r, 9, <j>) coordinate system is 
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7,r' 



x-l 



x-l 



with the poloidal component 
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2 lx-1 



2(x-l) 2 -x-l V*" 1 



-A6»b cos 
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(A4) 



(A5) 



where A9b = Bg /B r is the angle that a bending field line makes with the local radial vector. To evaluate the lab-frame radiation 
force ( l20| > we relate the solid angle of incoming photons to the emission coordinates via 



jdQ-ydcfy-y . 



(x- iy 

The poloidal and toroidal radiation force is then evaluated as follows. We begin by writing 

l-(3-k = A + B cos 4> 1 +C sin 4> 1 
and express the components of the wave vector as 

j3 Pt< j ) -k = Dp i< j > + E Pt< f, cos 7 +F Pi< j> sin 7 . 
Integrating first over 7 and then 9 1 at the 'emission' surface gives 



c{x-\y Jo 
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where G p and H p are accurate to first order in T" 2 and (9 jx) 1 . The equilibrium Lorentz factor of the photon field, the frame in 
which F md vanishes, is found by solving T' x = in (J47j). The results are shown in Figure [l5]for the poloidal field line profile 
described in Sec. 13.21 At a radius x-l 3> Oj one finds r eq ~ x/9j, with a coefficient of order unity that depends on the footprint 
angle and flaring profile. A thin jet defines a relativistic frame at relatively small distances from the 'emission' surface, as 
compared with a spherically symmetric outflow (for which r eq ~ 3 1,/4 x). Photons arriving at a scattering point from large angles 
provide relatively strong drag. 

The jet fluid maintains rapid rotation around the light cylinder, Of ~ 1 /xlu, where the fluid flow is less aligned with the radiation 
field and r eq is reduced. Estimating B^c^l /xuOf (just outside the light cylinder), one finds 



1 



3(X- l) 4 x4 w 404 
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f. 



valid for all x. 
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FIG. 15. — Lo rentz factor at which material in a thin jet (9j = 0.2) feels a vanishing radiation force F rad . Upper panel: Moderately flared magnetic field [59 = 1 
in equation 451H . with black lines corresponding to different footprints at the breakout radius: 8f t = 0.05,0.1,0.15 (solid, dotted and dashed). Lower panel: 
Different degrees of magnetic flaring, 56 = 10, 1,0.3 (weak, medium, strong flaring) for a field line anchored at 9f t = 0.1. Long-dashed red line: T e q for a 
spherical emission surface and monopolar poloidal field. 



ACCOUNTING FOR ROTATION OF THE PHOTON FIELD 



The photon source rotates rapidly in some cases, e.g. a rapidly rotating star such as a millisecond magnetar, or the merged 
remnant of a white dwarf binary. We can approximate the effect of a rotating emission surface by setting 



x6 



f 



(Bl) 



in equations d43l >. Here p^jt is a constant representing the aberration of the outflowing photons at r = r s (x= 1). In this situation, 
plasma near the emission surface can more easily co-rotate with the radiation field while still being accelerated outward. 

The value of (3^$ depends on the type of source. One has /3^ « ~ Qr s /c = lu when the photons flow from the surface of a star 
of radius r s through a transparent wind. On the other hand, if the outflow is optically thick in a narrow radial zone close to the 
engine, then one expects (3$^ ~ (f2r,/c) _1 ~ u~ l based on the conservation of angular momentum from the light cylinder out to 
the transparency surface (x = 1). 

WIND EQUATIONS FOR JET MODEL II 

Here we derive the equations (l46l-(l50l for the longitudinal development of Lorentz factor and particle angular momentum 
along magnetic flux surfaces. Beginning with the poloidal and toroidal components of the Euler equation, (BTt and (l42t . we 
expand the derivatives on the right hand side as 



d, (r6B^=B 4> d l {rd f ) + r6^diB r + r e f B r d, 

This can be evaluated using di(r9f) ~ 9f + Bg/B r = Oj + AOb, equation (f39l > for diB r , and 
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where 
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The wind equations (HTt , (l42l are now transformed into ordinary differential equations by ignoring the cross-field force balance 



B r V,T 3 

A6 B \ d,T 



n f (9 f+ A8 B ), 



(CI) 



(C2) 



(C3) 



1- 



(x6 f u;) 2 T 



1 B d 





dT 


i 


/ B r 




"^7 . 



(x6 f uj) 2 r By 



A 



dL m (T* Ba 

— = 2~r + x*Rj\ 



dl (x9 f uj) 2 B r 



(C4) 



Hot Electromagnetic Outflows II: Jet Breakout 



23 



(x9 f uj) 2 T 
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(r 3 >a). 



(C8) 



Two simple tests of these equations are made possible by neglecting the radiation force. The energy and angular momentum 
integrals (|38T > are now related by 

— Q (C7) 

dl dl 

This equation is recovered by summing (IC4) and B^jB r times dC5b . and making use of Ferraro's law (3). Second, outside the 
fast point the inertia of the magnetofluid is dominated by the matter: the coefficient of dT / dl in equation ( lC41 i is ~ 1 - <r/T 3 and 
approaches unity. Since in addition — > 0, the term involving C m can be neglected and one finds 

dT a* dA 
~dl ~ A 2 'dl 

The same result can be obtained from the integral equations (JHJ, ( l38T l and ( l39l >, 

r = r*+cr»(i-A _1 ). (C9) 

Equation d46b is then obtained by solving (IC4b and JC5t for dT /dl and dC m /dl. 

CROSS-FIELD FORCES 

Though we ignore the cross-field force balance, it is useful to estimate the transverse radiation force and compare it with the 
Lorentz force that is implied by a given field-line profile. Our procedure becomes inconsistent if the transverse radiation force 
dominates, because the radiation field will then comb out the field lines in the radial direction. 

The polar component of equation < f20b gives the estimate, 
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The first term on the right-hand side represents the force imparted by photons streaming from a finite polar cap toward particles 
on off-axis field lines. The second represents the drag imparted as the poloidal particle flow bends across the radiation field. The 
cross-field Lorentz force is given by 

mB 2 , „ a 



— [(V x B) x B] 6 

4irp 



AnprOf rdf 

Requiring this to be greater than (1DH gives an upper bound on the radiation compactness at jet breakout (x > x* 3> 1), 

.2 //i a/. \ -1 



X* < 2xT 
a x*6f 



Or A9 h 



(D2) 
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Our calculations can admit values of x* as large as ~ 10 4 a without any inconsistency, given the typical jet parameters 9j ~ 0.2, 

x*~3,r>r*~io. 
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